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SUMMARY 

An approximate method Is presented for calculating Idle pressure 
distribution on conical "bodies of noncircular cross section in a super- 
sonic flow field. By a superposition of elementary conical flows due 
to line sources, the flow about an arbitrary cone may be described. 
Illustrations of the pressure distribution about several shapes are 
included to demonstrate the method. The problem of such a body at 
angle of attack may also be solved by the same method, as well as the 
problem of yawed flight at angle of attack. 


INTRODUCTION 

Linearized methods have been applied to determine solutions of 
the supersonic flow field about both solid and open-nosed bodies of 
revolution. (See references 1 to 4.) The methods of reference 1 have 
been extended herein to produce a linearized solution for cones of arbi- 
trary cross section, such as might serve as forebodies of nonsymmetrical 
fuselages. The solution is based on the use of a combination of line 
sources inclined arbitrarily to the flow direction. If the proper, sources 
are chosen, any conical body shape can be described and the resulting 
surface pressures can be calculated. 

The method presented herein was devised during the fall of 1947 
at the NACA Cleveland laboratory. 


SYMBOLS 

The following symbols are used in this analysis : 
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constant proportional to source strength 
free -stream Mach number 

slope of line source with respect to flow direction 

component, normal to flow direction, of distance from point 
on line source to arbitrary point in flow field 

free-stream velocity 

component of total velocity normal to body at surface 

radial, (cylindrical coordinate) perturbation-velocity component 

axial perturbation-velocity component 

tangential (cylindrical coordinate) perturbation- velocity component 
cylindrical coordinates 
cotangent of Mach angle, Vm 2 - 1 

angular position of line source, measured from 6 - tt/2 plane" 

angle between radial velocity and normal to body, measured in 
x = constant plane 

angle whose tangent is ratio of radial to axial coordinates 
of point on body surface, tan"^- r/x 

axial coordinate of .line source 

source strength per unit of axial length 

perturbation-velocity potential 

angle between normal to body surface and normal to ray, 
which intersects axis of body 


GENERAL ANALYSIS 

In the solution of the pressure distribution over a body, a method 
of successive approximation is used. A perturbation-velocity potential 
based on linearized flow is found from which the three velocity components 
may be determined. If these velocities satisfy the boundary conditions 
for the desired body, the potential describes the flow about the actual 
body. The pressure coefficient may then be found from 
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C p = - 2U X /CJ 


( 1 ) 


The perturbation -Telocity potential may not be assumed com- 
pletely arbitrarily, but is subject to several general limitations 
and to some particular ones imposed by the body. The particular 
l imi tations form a guide in the selection of the potential. The 
two general limitations are: (1) the potential must satisfy the 

Prandtl-G-lauert equation; and (2) under the conditions of linear- 
ized flow, the disturbance velocities should disappear at the Mach 
cone. 

The Prandtl-G-lauert equation for compressible frictionless 
potential flow, in linearized font and cylindrical coordinates is 


(M 2 - 1) 


B^p l 5 
Bx 2 " r 



i qe - o 
Be 2 


( 2 ) 


A line source OT of strength per unit length f (|), lying 
in the plane 6 = [(it/2) + 6], whose slope with respect to the 
flow direction is m, is shown in figure 1. The potential at same 
point P due to such a disturbance is 


'x-R0 


cp = i 


f (l)d| 


*/ (| - x) 2 - p 2 B 2 


where R is, from the geometry of figure 1, defined by 
R 2 = r 2 + (m|) 2 - 2m£r sin (6 - 8) 


This potential can readily be shown by direct substitution to 
satisfy equation (2). 

Combining equations (3) and (4) and rearranging the terms 
result in 


(3) 


( 4 ) 




Die limits of integration are the origin 1=0, because there is no disturbance 
ahead of it and no disturbances are propagated upstream; and \ = x - HP, because no part 
of the line source beyond that point includes the point (x,r,9) within its Mach cone. 


Equation (5) is of the same form as equation (3), an d can be readily integrated if 
f(£) is assumed proportional to £ : that is, f(|) = ( - irrE^l - m 2 ??){ . Integration under 



The perturbation- velocity components U^, U x , and Ug, respectively, can be found by dif- 
ferentiating equation (6) with respect to r, x, and 0; 
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0> 


cos (6-5) 

— — - cosh*"-*- < 

[ 1 / sin (0-8)] 


- 2-2 
1-m P 

1j j [l- 5 ^ 2 sin (0-8)] - ^^(isaV) 


[l - ~^ 2 sin (0-5)] 


sin (0-5)] ? - (l-2-g!) (l-m 2 P 2 ) J 


„ 2 ft 2 
L-ii- 


fl-mV 


(9) 


If m = 0, that is, if the l ine source is in the flow direction, these equations reduce 
to thoBe developed by van Birman and Moore for the flow about a right circular cone (see 
reference 1). 



(7a) 

U x - -E cosh" 1 fp 

(8a) 

U a -o 

(9a) 


The boundary condition for a particular flow requires that the velocity normal to the 
body at the surface be zero. In figure 2(a), OTP is normal to the x-axls and Intersects It; 
0' G is normal to a ray 00* of the body and intersects the x-axis; 0‘H is normal to the 
body surface; O’L is normal to the contour of the body section found by cutting the body 
with a plane perpendicular to the x-axis and lies in that plane; and O’D is normal to plane 
OO'GEPft. 


965 


NAGA TN No. 1659 



MCA TN No. 1659 


7 


From figure 2(a), 


U = 0 = 
n 



cos u - (U + U ) sin v 

• 


cos - Ug 


sin 


* 


or 



tan \1/ / » 

% - - (U + U x ) tan v 

(10) 

But, 

also from figure 2(a) , 



tan 'k = tan X cos V 

(11) 

and. 

with reference to figure 2(b), 



tan X = 

rd0 

(12) 

where r = f (0) is the equation of the body cross section. 

If equa- 


tions (10) to (12) are combined, the boundary condition for the flow 
is obtained: 


°r - ? § *9 ' + ^ ! <“> 

The second general condition for linearized conical flow, stated 
previously, is that the perturbation velocities approach zero as the 
Ifech cone is approached, or as rp/x approaches unity. This con- 
dition is exactly true for equations (7) to (9). 

The potential found from the single-line source that has been 
considered until now is insufficient to calculate the flow about an 
arbitrary body. Because the potential due to one source satisfies 
the general limitations of the problem, however, a series of poten- 
tials due to a number of sources of various strengths and positions 
can be so added together that the resulting flow satisfies the bound- 
ary conditions for the particular body in question and thus can be 
assumed to be the desired flow. 

Four parameters are considered in selecting the source pattern: 
the ang ular position 5, the slope of each source relative to the 
flow direction m, the number of sources , and the strength of each. 


APPLICATION OF METHOD 

Because the perturbation velocities are functions of only r/x 
and 0, the resulting body contours are conical and only one body 
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section that is in a plane normal to the free -stream-flow direction 
need "be considered. The problem is thereby reduced to two dimensions 
and the bbdy is completely described by one plane, as in figure 3. 

The line sources thus appear in such a* figure simply as points. 

With the use of several sources of varying strength and position, 
the boundary condition, defined by equation (13), can be made to yield 
the flow about the desired body shape. Although such a solution is 
found largely by trial and error, some general rules for the selection 
of sources yielding a specified body shape can be established. The 
four parameters previously mentioned, that is, the angular posi- 
tion 6, the slope of the line sources m, the number of sources, and 
the source strength, must be kept in mind. 

1. The axes of symmetry of the body section should be noted, 
inasmuch as the sources must be symmetrically arranged relative to 
the same axes. 

2. The source nearest to a peak in the section should be nearer 
to that point than to any other point on the section because the per- 
turbation velocities due to the body are a maximum at the peak and, 
inasmuch as velocities due to the source increase as the distance 
from the source decreases, the point closest to the source has the 
highest velocity. The distance from a point on the body to the 
nearest source must therefore be equal to or less than the radius 

of curvature of the body section at that point. 

5. If the body is elongated, a series of sources in a line are 
required. The slenderer • the section, the more sources are needed 
to prevent contour irregularities in the body described by the 
re axil ting solution (equation (13)). Also, for a given body of this 
type, as the section narrows the sources should be closer together. 

4. In general, the larger the number of sources used, the more 
closely the linearized flow- obtained fits the body in question and 
the fewer 'are the number of trial solutions required to obtain a 
satisfactory answer. Each such solution, however, is more laborious 
than one using fewer sources. 

To summarize, the angular position of the sources Is determined 
by rules 1 and 2, the distance from the axis by 1, 2> and 3, and the 
number of sources by 1 and 3. Rule 4 serves as an. over-all guide. 

Only E, which Is proportional to the source strengths, remains to 
be determined. From the symmetry condition, sources that are In the 
same position relative to the axes of symmetry -have the same strength. 
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With a known number of different source strengths, equation (13) may- 
be solved for the strengths at the s ame number of points on the surface 
by using equations (7) to (9). The potential is now completely defined. 
This potential should then be checked in equation (13) at several points 
on the body to determine whether the flow due to the - source potential 
is the same as that over the actual body. It should be noted that 
the strength of a source may be negative. 

The case of flow at angle of attack can be solved by considering, 
at zero angle, a body whose cross section consists of sections of the 
actual body taken normal to the free stream instead of normal to the 
axis. The x-axis of the new body is then par all el to the flow direc- 
tion. For a small angle of attack, the two cross sections probably 
differ little, but their positions relative to the coordinate axes 
differ. The case of yawed flight at angle of attack can be solved 
in a similar manner (fig. 4). 


Examples 

Several examples follow to illustrate the general rules that 
have been outlined: 

Example I. - Assume that an elliptic section (fig. 3) is desired. 
Frcm the symmetry condition, two sources are assumed as a first approx- 
imation to such a section. These sources are of equal strength, equi- 
distant from the x-axis, and at 5=0 and 6 = jt. The velocities 
are found from equations (7) to (9). 

The equation of an ellipse is 

2 d 2 c 2 

r ~ 2 2 2 2 
d sin 6 + c cos 6 

where c and d are major and minor semi -axes, respectively. 

1 dr 

From this relation, — g-g may be obtained: 


1 dr _ sin 6 cos 9'(c^-d^) 


r d© 


d sin 9 + c cos 9 


(14) 


When a Ifech number and values for c and d, are given, a position 
for the sources (that is, a value of m) is selected. From, the 
second general rule, c - m of figure 3 should be somewhat less than 
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•the radius of curvature at c, «/2. With this guide, a value of m 
can he assumed. From the coordinates of one point on the surface K/U 
can he found by substitution of equations (7) to (9) and (14) in equa- 
tion (13). Several other points should then he checked in the boundary 
relation. If agreement is poor, the sources should be moved and per- 
haps others added. When a potential is obtained that satisfies equa- 
tion (13) closely enough, the pressure coefficient can then be found 
from equation (l) . 

An example of the results of such a calculation is shown in 
figure 3. The body shape obtained by using the value of K previ- 
ously determined and by solving equation (13) for r/x as a function 
of 0 is shown together with the desired shape. In this example, where 
the contour is almost circular, the pressure distribution approximately 
follows the body shape. The deviation of the calculated section from 
the desired ellipse could be considerably decreased by placing the two 
sources farther apart and adding a third source at the origin. 

If the flow about an ellipse having a larger ratio of major to 
minor axes than was used In figure 3 is desired, the perturbation- velocity 
potential for two sources gives poor results. The effect of varying K 
while holding m constant is shown in figure 5. A similar result is 
obtained by holding the length of one of the axes constant and varying 
the length of the other by changing m, while holding the source strength 
constant. In order to obtain a satisfactory solution for an ellipse having 
the ratio of the axes much greater than that in figure 3, a series of sources 
can be used. 

Example II . - Assume that an ellipse having the ratio of major to 
minor axes equal to 3 is desired. (See fig. 6.) Obviously, for symmetry, 
a number of source pairs, as used in the previous example, plus perhaps 
a single source at the origin, will give the desired solution. The position 
of source 1 in figure 6 may be assumed from the condition that its distance 
Xrcm the peak must be approximately the radius of curvature at the peak. 

The distance between sources 1 and 2 is taken as about equal to the distance 
from source 2 to the nearest point on the body. The remaining sources may 
be similarly chosen. This procedure gives a system of three source pairs 
whose velocities are found from equations (7) to (9), and a single source 
whose velocities are found from equations (7a) and (8a), which are simply 
the velocities found when m is zero. When equation (13) Is solved for 
this system of sources at four points on the body, K, , M, Kg, and K4 
(fig- 6) are determined and the pressure coefficient can be found from 
equation (l). 

The desired ellipse a n d the contour calculated by the use of seven 
sources are shown in figure 6 together with the pressure distribution 
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corresponding to the source configuration. In this example the section 
deviates greatly from a circle and the pressure distribution no longer 
follows the body contour. The maximum and minimum values of the pres- 
sure coefficient, however, still occur at the maximum and minimum points, 
respectively, on the section. 

Example III. - Now assume that the pressure distribution is desired 
over a body whose cross section is a triangle modified by rounding the 
vertices (fig. 7). From symmetry considerations, three sources of equal 
strength should be assumed at 5=0, 2it/3, and 4it/3. These sources 

should be placed at approximately the center of curvature of the vertices. 
Because the body is approximately circular, a source that is not of the 
same strength as the others should be placed at the origin. 

The velocity components can then be calculated from equations (7) 
to (9) and the strength of the sources can be established by solving 
equation (13) at two points. 

Such a surface, with the corresponding pressure distribution, is 
illustrated in figure 7. As in the first example, the pressure distri- 
bution follows the trends of the body shape. 

It must be remembered that the examples given are meant to illustrate 
the method of solution of such bodies rather than to show actual pressure 
distributions, although the trends indicated should be correct. The bodies 
chosen are probably not slender enough for great accuracy in a linearized 
solution. 


SUMMARY CF ANALYSIS 

An approximate method has been presented for calculating the 
pressure distribution on conical bodies of arbitrary cross section in 
supersonic flow. By a combination of elementary conical flows due to 
line sources, the flow about a slender arbitrary cone can be described. 
Four parameters are considered in determining such a system of soirees : 
the spacing of the sources around an axis lying in the flow direction, 
the slope relative to the flow, the number of sources, and their strength. 
The first three parameters can be determined by several conditions. 

First, the same symmetries will hold for the sources as for the body. 
Second, the distance from a peak to a source will be less than or equal 
to the radius of curvature at that point. Third, for an elongated body, 
the slenderer the section, the closer together the sources must be. 
Finally, the greater the number of sources used, the more accurately 
the desired body can be approximated. The solution will be more labo- 
rious, however, if more sources are used. The fourth parameter, source 
strength, may be found by direct calculation. 
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The case of a conical "body at angle of attack and of yawed flight 
can he solved by the same method. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 27, 1948. 
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■- Desired section 

Calculated section 

Pressure distribution 

O Source position 



Figure 5. — Comparison of desired elliptical body with calculated body and 
resulting pressure distribution. (3 = 1.5; m = 0.125; K/TT “ 0.0220. 
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Figure 5. - Surfaces obtained with sources of varying strength 

at (a, ± n/ 2 )» m = 0.20. 
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“ 0.26 

Kj/0 - 0.00483 


■s - .i» 

Ka/U - .00396 

— — O.loulat.d body 

«3 - .10 

Kj/tT * .00361 

Pr.is'or. dl.trlbution 

■ 4 - 0 

Kj/tT « .00327 

O Sours, position 



Figure 6* — Comparison of desired alender elliptical body with calculated body and resulting 

pressure distribution, p ^ 1.5. 
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— — — ■ Desired and calculated body 

1 - -——Pressure distribution 
O Source position 



Figure 7. — Comparison of desired modified triangular body with calculated 
body and resulting pressure distribution, p «= 1,5; m = 0.16; = 0.0102; 

ICg/tr = 0.0154. 




